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Abstract 

We re-examine a recent proposal for obtaining f3 using the measurements of the 
b —> d penguin decays B%(t) — > K°K° and B®(t) — > (j)K s , along with a theoretical 
assumption. We show that there are in fact three assumptions one can make, so 
that the method can in principle be used to extract a, f3 or 7. We also show that it 
is the assumption which yields 7 which is the best. However, the theoretical error 
on this assumption is still 25-30%, which leads to an error on 7 of at least 20-25%. 
Given our current understanding of hadronic physics, it does not seem possible to 
reduce this error. 
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Measurements of CP-violating rate asymmetries in the neutral B system will 
allow one to obtain the CP angles a, (3 and 7 [![]. From these the unitarity triangle 
can be constructed, and one can test the predictions of the standard model (SM). 
If we are lucky, these measurements will reveal the presence of physics beyond the 
SM. 

One of the best ways of searching for new physics is to consider two different 
decay modes whose CP asymmetries probe the same CP phase within the SM. Any 
discrepancy between the measured values of these asymmetries will point unequiv- 
ocally to the presence of new physics. One possibility is to consider pure b — > d 
penguin decays such as B d (t) — * K°K° or B®(t) — * 4>K S . At the quark level, these 
decay amplitudes take the form b — > dss. If such decays are dominated by inter- 
nal t-quark exchange, the amplitude is proportional to V t bV t * d , where the are 
elements of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix. In the 
(approximate) Wolfenstein parametrization of the CKM matrix ||, the only two 
matrix elements which have a nonzero weak phase are V id [~ exp(— if3)\ and V u b 
[~ exp(— 27)]. Thus, CP asymmetries in pure b — > d penguin decays probe the weak 
phase (3. By comparing this value of (3 with that extracted via the conventional 
mode B d (t) — > J/ipK s , one can search for new physics. 

Unfortunately, the b — > d penguin amplitude is not dominated by t-quark ex- 
change. For example, the uu quark pair of the tree-level decay b — ► uud can rescatter 
strongly into an ss quark pair, giving an effective V u bV* d contribution to the b — ► d 
penguin amplitude, and similarly for the b — > ccd tree-level decay. Buras and Fleis- 
cher |3J] have noted that the u- and c-quark contributions can be between 20% and 
50% of the leading t-quark contribution to the b — > d penguin amplitude. And since 
V u bV* d and V c bV* d have different weak phases as compared to V t bV* d) this implies that 
CP asymmetries in pure b — > d penguin decays do not cleanly probe the weak phase 
(3. 

But this then begs the question: is it possible to isolate the t-quark contribution 
to the b — > d penguin amplitude? If so, one could then obtain (3 from this piece of the 
amplitude, and compare it with the value found in B d (t) — > J/ipK s . Unfortunately, 
as was shown in Ref. [||], this is not possible. In a nutshell, the argument is as 
follows. The b — > d penguin amplitude can be written generally as 

p = PuV ub v: d + p c v cb v: d + p t v tb v; d . (1) 

Now, the unitarity of the CKM matrix implies that V ub V* d + V cb V* d + V t bV td = 0. 
Thus, any of the three pieces in the above amplitude can always be eliminated in 
terms of the remaining two. For example, if we eliminate the V u bV* d piece, we obtain 

p = (p c -p u )v cb v: d +(p t -p u )v tb v td 

= V cu e l& ™ +P tu e i5 ^e^ , (2) 

where we have explicitly separated out the weak and strong phases and absorbed 
the magnitudes |V^,T4d| an d |V^Vtd| into the definitions of V cu and Vtu, respectively. 
Similarly, eliminating VtbV t * d gives 

P = V ct e iSct + V ut e iSut e" i7 . (3) 
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Now, suppose that a method existed which would permit one to extract the CP 
phase (3, with no hadronic uncertainties, using the parametrization of Eq. (0). If so, 
then (3 could be expressed entirely in terms of experimentally measured quantities. 
However, Eq. (|3[) has the same form as Eq. (fj). Thus, this same method would allow 
us to cleanly obtain —7 using Eq. (|). In particular, —7 would be expressed as the 
same function of observables as /3, leading to the conclusion that f3 = —7. Since this 
is clearly not true in general, one concludes that it is not possible to cleanly measure 
the weak phase of the i-quark piece (or indeed any other piece) of the b — ► d penguin 
amplitude. In Ref. 0], this is referred to as the "CKM ambiguity." 

However, in Ref. || it is also argued that it is possible to resolve the CKM 
ambiguity, and hence isolate a particular piece of the b — > d penguin amplitude, if one 
makes an assumption regarding the hadronic parameters involved in this amplitude. 
This fact has been used by two of the present authors (Kim, London) and Yoshikawa 
in Ref. || to obtain (3 from the t-quark piece of the b — > d penguin amplitude. 
The idea is to use the two penguin decays B d (t) — > K°K° and B®(t) — > (j)K s . If 
one eliminates the V u bV* d piece from the penguin amplitudes, as in Eq. (|), the 
amplitudes for these two decays can be written as 

A(B° d -> K°W) = V cu e lS ™ + V tu e l5t - e~ ip , 
A(B° s ^<f>K s ) = V' cu e lS ™ + V' tu e i5 - e~ ifi . (4) 



The following assumption is now made: 

'"■u 'Pcu/'Ptu 



(5) 



tu 



(Note that the dependence on the CKM matrix elements cancels in this ratio, so that 
this really is an assumption about the hadronic parameters of the two amplitudes.) 
Measurements of B d (t) — > K°K° and B®(t) — > 4>K S , combined with the assumption 
in Eq. (Bj), will then allow the extraction of f3, the weak phase of the t-quark piece 
of the b — > d penguin amplitude. By comparing this value of f3 with that found in 
B d (t) — > J/ipK s , one may be able to detect the presence of new physics. In what 
follows, we will refer to this as the "KLY method." 

However, there is a problem with this method which has been overlooked in 
Ref. ||. It is similar to the CKM ambiguity: how do we know that V u bV* d is the 
correct term to eliminate? For example, had we eliminated the V t bV t * d piece, as in 
Eq. (H), the amplitudes would take the form 

A(B d -> K Q K°) = V ct e i5ct + V ut e i5 ^ e* 7 , 
Aift^+Ks) = V' ct e iS '« + V' ut e<<e^ . (6) 



Using the assumption that 



n _ Pa/Pi 



ut 



(7) 

r't V'jV'ut ' ' 1 ) 

the KLY method can be used to allow us to obtain 7. Similarly, the elimination of 
the V c bV* d piece in the amplitudes, combined with the assumption that 

_ 'Puc/'Ptc _ , /q\ 
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would permit the extraction of a. The confusion regarding the validity of the three 
assumptions in Eqs. (|5|), (0) and (H) can be thought of as a "second CKM ambiguity." 

It is clear that the three assumptions in Eqs. (§), (0) and (§) cannot be simul- 
taneously true. If they were, then, as was the case in the discussion of the CKM 
ambiguity, one would obtain results such as (3 = —7, which do not hold in general. 
We therefore deduce that (at least) two of the assumptions are poor. That is, the 
values of the CP phases obtained using these assumptions will differ enormously 
from their true values. 

In fact, it can be argued that all three assumptions are likely to be poor. Af- 
ter all, even though they are based on the same quark-level process, the hadronic 
quantities (form factors, strong phases, etc.) describing the decays B®{t) — > K°K° 
and B®(t) — > (pK s should be quite different. Thus, there is no reason to expect any 
of the ratios in Eqs. (||), (^) and (f|) to equal one, in which case it appears that the 
method is of little practical use. 

However, this is not as serious a problem as appears at first glance. Consider 
again the assumption of Eq. ([|). In fact, the KLY method does not require that 
r u /r' u = 1. All that is necessary is that we know the value of this ratio. Thus, if we 
could theoretically calculate the value of r u /r' u , the KLY method could then be used 
to obtain (3. Similarly, if we could calculate r t /r' t or r c /r' c , we could extract 7 or a. 
We therefore see that the second CKM ambiguity is not necessarily a disadvantage: 
depending on what hadronic information is known, any of the three CP phases could, 
in principle, be obtained using the KLY method. 

This is therefore the key question: how well can we estimate the values of the 
three ratios r u /r' u , r t /r' t and r c /r' c 7 Can we make general statements regarding this 
question, or is it process-dependent? Also, for a given theoretical uncertainty on a 
particular ratio, what is error on the corresponding extracted CP phase? These are 
the issues which we address in this paper. 

We argue below that, with our present theoretical understanding of hadronic 
B decays, neither r u /r' u nor r c /r' c can be computed with any degree of reliability. 
Therefore the KLY method cannot be used to obtain j3 or a. On the other hand, for 
pure b — > d penguin decays such as B®(t) — > K°K° and B®(t) — > (pK s , r t /r' t can be 
calculated to lie in a much narrower range. Thus, for these decays, the KLY method 
can be used in principle to extract 7. Unfortunately, as we will see, the error on 7 
remains fairly large, in the range of 20-25%, so that this method cannot be used to 
obtain a precise measurement of 7. We examine possible ways to reduce this error. 

We begin our analysis by addressing the question of how to estimate the ratios 
r u /r' u , r t /r' t and r c /r' c . The important point to realize is that the Pj's which appear 
in the ratios [Eqs. @, (^), (H)] are actually matrix elements of penguin operators. 
Thus, in order to answer this question, we need some sort of framework in which to 
evaluate hadronic matrix elements. Factorization is usually employed to calculate 
nonleptonic amplitudes. Corrections to naive factorization have been calculated 
within QCD factorization [7|] and the perturbative QCD (pQCD) approach. In 
QCD factorization, naive factorization is recovered as the leading-order term, and 
one systematically computes corrections to it in an expansion in a s (mt>) ~ 0.2 and 
AQco/^b- At 0(a s ) the vertex and hard spectator corrections modify only the top 
penguin amplitude and do not introduce any additional weak phases. The penguin 
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or the rescattering corrections generate the up and the charm penguin pieces. The 
SM effective hamiltonian for hadronic B decays is 0: 

Hlff = ^[v fb v; q ( Cl oi f + c 2 o« f ) 

10 

- EW* W + V cb V^ + VaV^O?} + h.c. (9) 

i=3 

Here, q can be either a d or an s quark, depending on whether the decay is a AS = 
or a AS = — 1 process. In the first terms, / can be a u or a c quark, while in the 
last terms, the superscript u, c or t indicates the flavour of the internal quark. Note 
that, in the effective Hamiltonian H e ff above, we have explicitly included the up 
and charm penguin pieces which are generated by rescattering. The values of the 
Wilson coefficients c{ for the penguin operators evaluated at the scale jj, — rrib — 5 



GeV, for m t = 176 GeV and aJmz) = 0.117, are 10 



c* = 0.017 , 4 = -0-037 , 4 = 0.010 , 4 = -0.045 



\ = -1.24 x 10" 5 , 4 = 3 - 77 x 10 " 4 > 4 = -0.010 , 4o = 2- 06 x 10 ~ 3 > 

4,5 = -c\JN c = p:/n c , 4 9 = pi , 4 10 = , i = u, c , (10) 



where iV c is the number of colors. The leading contributions to P % se are given by 
PI = (t) c 2(f + %M 2 )) and = (^)(iV cCl + c 2 )(f + GK, g 2 )), where 
the function G(m, fi, q 2 ) takes the form 



G(m, /x, g ) = 4 / x(l-x)ln ^ — dx , (11) 

JO /i 2 



where q is the momentum carried by the virtual gluon in the penguin diagram. 

Of course, as mentioned above, we are really interested in the matrix elements of 
the various operators for the decay B — > fif 2 - We therefore define new coefficients 



—u,c 

c 4 as 





cT C (Q 2 )Oi \B) 


(hf2 


Oi\B) 



-u,c = wi V, ^ p (12) 



The values of c"' c can be calculated in the approaches of Ref. and Ref. f| if 
the light cone distributions of the various mesons are known. However, the values 
for the 4' C wm ! m general, be different in the QCD factorization and the pQCD 
calculations. This is because pQCD assumes that the the light quarks forming the 
final state light mesons must all be energetic while QCD factorization assumes that 
the spectator quark coming from the B mesons remains soft as it combines with 
an energetic quark to form one of the final-state light mesons. Without adopting a 
particular approach, we will follow the usual practice, which is to simply replace 

-u,c _^ c «> c ( g 2j ^ ( 13 ) 

and we will conservatively allow q 2 v to vary between m 2 /4 — > m 2 /2. (Note that this 
is the range which has been used in the past to take into account possible process 
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dependence |12|.) With this prescription, the relations between the various c 
are the same as those between the various c"' c (g 2 ). In particular, we still have flTB" 



— u,c — u,c 



(14) 



In order to calculate matrix elements, we need to consider specific final states. 
Consider first the decay B® — > K°K°. Using the naive factorization approximation, 
along with the fact that Cg' c = [Eq. (|i4|)1, one can write 



# c (l - W ) [(Oll) ~ 2 (0 SP )] 



(15) 



where 



(Oll) 
(Osp) 



K° 
K° 



S7m(! -7s)6 



s(l- l5 )b 



B° d )(K° 



B°) K 



dy(i- 75 )s|o) 

d(l+ 7 5)s|0) , 



(16) 



and we have dropped factors common to P u , c ,t- (The operator Osp appears due to a 
Fierz transformation: (V-A)® (V + A) = -2(S - P) ® (S + P).) The contribution 
from the top penguin amplitude is given by 



Pi 



2 lCl0+ N r 



(Oll) - 2(4 + (Osp) 



In the above, we have neglected the contributions from cj-, 
It is convenient to rewrite P u c and Pt as 



(17) 



where 



and 



- j^)X r (Oll) 



P 



+(a 4 - a 6 - -aio) (Oll) + a%X r (0 LL ) 



Ci+°p , 2 = 4,6,10 , 
Q + ^ , i = 3, 5, 9 , 



(18) 
(19) 
(20) 

(21) 

m b m s + m d _ 

Note that, for m K = 500 MeV, rrib ~ 5 GeV, m s ~ 100 MeV and ~ 0, one finds 
that X r ~ 2 for this decay. 

Another possible decay mode is i?^ — ► K*K*. In this case the above analysis is 
unchanged, except that 



2 (Osp) 
(Oll) 



For B?, -> K ^ , 



A" r 



2m 2 

1 + 



AV 



(22) 
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since here (Osp) = 0. (Note that the KLY method does not apply to the decay 
— ► KK* since this final state is not CP self-conjugate.) Finally, one can also 
consider decays such as B® — > (j)K s . However, over a large region of parameter 
space, P u = P c ~ JCJ, so that the CP asymmetry probes /? directly, and the KLY 
method does not apply. We therefore concentrate only on B® — > K^*'K^*' decays in 
the analysis below. 

One can now construct the ratios r u , r c and r t defined in the numerators of 
Eqs. (|5|), (|7|) and (|§|). We use r u0 , r c0 and r t0 to denote the values of r u , r c and r t in 
the naive factorization limit. In addition, it is useful to separate the dependence of 
r uo, r co and r t0 on the CKM matrix elements from the hadronic dependence: 



v; b v cd 



V t tV td 



(had) 
uO 



V* b Vud 



v t tv td 



(had) 
cO 



no 



v; b v cd 



V ub V ud 



Shad) 

to 



(23) 



Using the above expressions for the Pi, we find 



(had) 


Pc 


-Pu 




x r {\ - 


N 2 C ) (c§ - 






uO = 


Pt 


-Pu 




a 4 + a 6 (X r — 1) 


— 2 a 10 — 


x r {\- 


Nl)c% 


(had) 


Pu 


-Pc 




X r (l ~ 


- Nl) 


-cl) 




cO — 


Pt 


-Pc 




(a 4 + a e (X r - 1) 


— \ a io) 


-x r {\ 


- met 


(had) 


Pc 


-Pt 




(a 4 + a 6 (X r - 1) 


— 2°lo) 


-x r (i 


- Nin 


to — 


Pu 


-Pt 




(a 4 + a 6 (X r - 1) 


— 2°lo) 


-x r (i 


- 



(24) 



One obtains similar expressions for the ratios r 



/(had) 
uO ' 



j(had) 



cO 



i /(had) 

and r t0 



where the 

r- (i = u,c,t) are defined in the denominators of Eqs. (H), (0) and @. Note that 

the various ratios r„o°d) to an d r uo?co]to depend only on the Wilson coefficients, on q% v 
and on the quark masses in the factor X r . There is no dependence on hadronic 
quantities such as form factors and decay constants since they cancel in the ratios. 

So far we have considered only the rescattering corrections which occur at a s , 
and which generate the up and charm penguins. However, as noted earlier, there are 
additional vertex and the hard spectator corrections, which are also 0(a s ). These 



can be taken into account by the replacement 



eff 



ai(l+ti) in Eq. flU), where 



ti ~ 0(a s ) are process-dependent corrections to the naive factorization assumption. 
Since the corrections ti depend on several poorly-known nonperturbative quantities, 



we will treat them as free parameters. The process dependence of the 



Jhad) 
u,c,t 



and 



then comes from three sources: (i) the value of the momentum transfer g 



2 



r(had) 
^"u,c,t 

in Eq. (|T3"D, (ii) ti, the 0(a s ) corrections to the aj, and (iii) the dependence of the 
quantity X r [Eq. (|20D l on the quark and hadron masses. 

Including now all the corrections to naive factorization, to first order in a s and to 
leading order in Agc^/mf,, we can make the replacement Pt — > Pt{l +x), where the 
process-dependent quantity x ~ a s (m b ) ~ 0.2 4 . We can then obtain the corrected 
values of r^ iad \ r^ hadS> and rf' ad ^ as 



„(had) 



(had) 
uO 



1 + ^OT + 2a; ^ COS ^-^ 



\Ptu\ 



4 We have neglected a possible small complex phase in x. 



6 



.(had) _ Jhad) 1 

c — cO / == 

^1 + x 2 ^ + 2xj^ cos(6 tc - St) 
1 + x 2 j^ + cos(5 te - 5 t ) 



.(had) r (had) + X \K¥ + Tfici CQS l° te ~ °t) 

rt0 i i+x2 w +2x \tu c ° s ( Stu - s J 



(25) 



where P t — P u = \Pt\e l5t — \P u \e l5u = \P tu \e lStu , and similarly for P t — P c . As usual, 
there are similar expressions for r'^ had ^ , r'j had ^ and r' t °^ • From these expressions one 
can see that the nonfactorizable corrections tend to cancel in r\ ; ^ and r't had \ and 
so these ratios are the least affected by such effects. Note also that, as mentioned 
previously, the dependence on the CKM matrix elements cancels in the ratios fj/r^, 
i.e. 



(had) 



r .i l(had) 
1 i T; 



(26) 



We are now in a position to calculate the theoretically-allowed ranges of the 
three ratios of Eqs. (|5[), (0) and (^j). There are several factors which can contribute 
to these ranges: since the quantities rf 10 ^ and rf had ' (i = u, c, t) are calculated for 
different processes, the momentum transfer q 2 v , the parameter X r [Eq. fl2TJ|)], and the 
nonfactorizable correction x may all be different. We therefore adopt the following 
procedure. For each of rf 1 ad ^ and rf had \ we allow the values of the quark masses 
to vary in the following ranges: 4.3 < m& < 4.9 GeV, 1.20 < m c < 1.30 GeV, and 
0.080 < m s < 0.120 GeV. We fix m d = 6 MeV. All masses are taken to be at the 
6-quark mass scale. In addition, for a given value of m&, we vary q\ v between m\j '4 
and ml/2. Finally, we take x to be in the range —0.2 < x < 0.2. 

In order to calculate the rf 1 ad ^ and rf had \ we must choose two decay pro- 
cesses. As a first example, we consider B° d (t) -> K* (892)^* (892) and Bg(t) -> 
tf*(1410)#*(1410). The quantum numbers (J PC ) of the if* (892) and if* (1410) 
are the same, and the latter can be interpreted as a radially excited K*. From the 
analysis of Ref. |T4| we can expect the branching ratio of B® — > i^*(1410)A"*(1410) 



to be comparable or enhanced relative to B® — > X* (892) X* (892) (]. (Note that since 
the final state consists of two vector mesons, an angular analysis will have to be 
performed to separate the three helicity states However, that does not affect 
our analysis here.) We find that r^ ad ^ lies in the range 0.14 < r^°^ < 0.54, and 

similarly for r'^ had \ Thus, we have 0.26 < r u /r' u < 3.86. Similarly, the range of 
r (had) and r i(had) ig from q 13 tQ 0>47) which leadg tQ q 28 < r j r ' c < 3.62. The al- 
lowed ranges for both r u /r' u and r c /r' c are clearly enormous. Since the KLY method 
requires a reasonably accurate theoretical prediction of these ratios, we therefore 
conclude that (3 and a cannot be obtained using this method. 

On the other hand, the allowed range for r^ had ' > and r[ } had ^ is considerably nar- 
rower: 0.92 < r\ had \ r't < 1.23. This is due partly to the fact that P u and P c 
are both quite a bit smaller than P t , so that the numerator and denominator of 



5 It may also be useful to consider the process B° d (t) -»■ if*(1680)AT*(1680), since the K*(16&0) 
has a significantly higher branching ratio to Kit (~ 39%) than the i^*(f4f0) (~ 7 %). This makes 
reconstructing the K*(1680)K *(1 680) final state easier than K*(U10)K*(U10). 
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r (had) roughly equal (and similarly for r' t ), and partly to the fact that the 
nonfactorizable corrections approximately cancel in r[ had ^ and r'} h [Eq. (p5|) ]. We 
therefore find that 0.75 <r t /r' t < 1.34, a considerably tighter range than that found 
for r u /r' u and r c /r' c . Thus, of the three ratios, r t /r' t has by far the narrowest range, 
so that in fact it is the CP phase 7 which can be extracted with the smallest error 
using the KLY method. Note also that the average value of r t /r' t in its range is 1.04, 
which is quite close to unity. Thus, the assumption of Eq. (|7|) is justified, though of 
course the key question is the error on the assumption. From now on, we therefore 
consider only the measurement of 7 using the assumption of Eq. ([/]). 

Another pair of processes that one can consider are B d (t) — > K°K° and B d (t) — > 
fr°(1460)K°(1460) (fT°(1460) has the same quantum numbers as the K°). Since 
these final states consist of two pseudoscalars, one must use the expression for X r 
found in Eq. fl2ip. However, the results do not change much: we find 0.94 < 
r[ had \ r' t [had) < 1.22, leading to 0.77 < r t /r' t < 1.30. 

Note that we have been extremely generous in estimating the allowed range or 
r t /r' t . We have allowed each of r[ had ^ and to take any value in their allowed 

ranges. That is, we have assumed that the momentum transfer q 2 av and the non- 
factorizable correction x can each take completely different values in the two decay 
processes. However, in practice, this is not likely to be the case if two similar final 
states are chosen. For example, we expect a similar value of q 2 av in the two de- 
cays B%t) -> #*(1410)K*(1410) and B° d (t) -»■ K°( 1460)7?° (1460), since the masses 
of the particles in the final states are almost equal. Also, it is reasonable to ex- 
pect that the effect of nonfactorizable contributions will be similar for the decays 
B° d (t) -> K* (892) K* (892) and B° d (t) -> K*(U10)K*(U10), since both final states 
consist of two vector mesons. Thus, it is likely that the range of r t /r' t may actually 
be much smaller than that calculated above, particularly for similar final states. We 
will come back to this point later. 

Of course, it is not enough to have established that it is 7 which can be extracted 
with the smallest error using the KLY method. What we really want to know is: 
what is the size of the theoretical error on 7 in this method? This is the question 
we now address. 

We first recall how measurements of two processes, B d (t) — > M X M 2 and B d (t) — > 
M[M 2 , along with an assumption about the ratio of penguin amplitudes, can be used 
to obtain 7. Using the convention B d = bd, we write the amplitude for B d — > M1M2 
using the parametrization of Eq. @: 

The amplitude for B° d -> M[M' 2 can be written similarly: 

Af^=V' ct ^+V' ut e^e^ . (28) 

The amplitudes for W d -> M X M 2 and W d -> M[M' 2) respectively Af lM2 and Af [M '\ 
can be obtained from the above amplitudes by changing the sign of the weak phase 
7- 

From time-dependent measurements of the process B d (t) — > M1M2, one can 
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obtain the following three observables: 



X 
Y 



aMiM 2 



+ 



M 2 



aMiM 2 



-Im f e -^ A MiM 2 * A M lM2 



+ Kt + ^ctV u t cos A cos 7 , 
VPetPut sin A sin 7 , 

= V 2 t sin 2(3 + 2V ct V ut cos A sin(2/3 + 7) 



+P 2 t sin(2/? + 2 7 ) 



(29) 



where A = 5 r t — 5,,+. One can also define a fourth observable: 



Z R = Re (e~^A^ M2 *A^ lM2 ) = V 2 t cos 2/3 + 2V ct V ut cos A cos(2/5 + 7) 

+P 2 t cos(2/? + 2 7 ) . (30) 

Given that the width difference between B® and B® is very small, it is unlikely that 
Z R can be measured experimentally However, note that Z R is not independent of 
the other three observables: 



Y' 



(31) 



Thus, one can obtain Z R from measurements of X, Y and Z r , up to a sign ambiguity. 
It is also useful to further define "rotated" observables: 

Z r = Zr cos 2/3 — Z R sin 2(3 

= V 2 ut sin 2 7 + 2V ct V ut cos A sin 7 , 



Z R = Z j sin 2(3 ~\~ Z R cos 2(3 

= V 2 ut cos 2 7 + Vl + 2V ct V ut cos A cos 7 . 



(32) 



Assuming that (3 is known independently (e.g. from the CP asymmetry in B — > 
J/ipK s ), we can obtain Z f and Z R from measurements of B®(t) —>■ M\Mi. The 
observables X', Y', etc. for the second process B® — > M^M^ can be defined similarly 
using the "primed" parameters of Eq. ( p8|) . 

Note that, apart from the CP phase (3, the three independent observables X, Y 
and Zj depend on four unknowns: V ut , Va, A, 7. Thus, one cannot obtain CP-phase 
information from the process B®(t) — > MiM 2 alone. However, the above equations 
can be solved to yield V ut and V ct as functions of 7: 



' ut 



X 



V. 



cos 27 — 1 

Z R cos 27 + Z l sin 27 — X 



ct 



cos 27 — 1 

Thus, combining both processes, one can write 



(33) 



n 



Vct/V„ 



ut 



Z R cos2 7 + Zj sin2 7 - X Z' R - X' 
Z' R cos 2 7 + Z'j sin 2 7 - X> ~Z~ h 



X 



(34) 
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O 1 /I o 
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27.0 


O O AO 

- 32.9 




a o 

0.2 


or /i o 

25.4 


- 3o. 3 




0.25 


24.4° 


- 38.3° 


P AO 

60 


A A 1 

0.01 


r a ^70 

59.7 


- 60.6 




a a r 

0.05 


r ^7 AO 

57.9 


- 62.4 




A 1 

0.1 


r r to 

55.7 


/I AO 

- 64.9 




A O 

0.2 


CT 1 C O 

51.0 


- 70.1 




0.25 


49.6° 


- 72.9° 


80° 


0.01 


79.6° 


- 80.7° 




0.05 


77.6° 


- 82.7° 




0.1 


75.0° 


- 85.4° 




0.2 


69.9° 


- 89.3° 




0.25 


67.4° 


- 93.5° 



Table 1: The range for 7 calculated from Eq. (^) assuming that r t /r' t = 1.0 ± Ar, 
for the input parameters of Eq. (^), and for three input values of 7 (■y in = 30°, 60°, 
80°). 

We therefore see that a prediction for r t /r' t will allow us to obtain 7. 

As shown earlier, r t /r' t is expected to lie in the range 0.75 < r t /r' t < 1.34. 
Although this range is far more narrow than those found for r u /r' u and r c /r' c , it is 
still very large. How does this translate into an error on the extracted value of 7? 
To examine this question, we take the true values of the theoretical parameters to 
be: 

Vet = 1.1 , Put = 1.0 , V' ct = 1.5 , V' ut = 1.36 , 

A = 30° , A' = 110° , p = 20° . (35) 

We also consider three values for 7: 30°, 60° and 80°. Given these inputs, we can 
calculate the values of the experimental quantities in Eq. (|34|). Then, given a value 
of rt/r f t , we can solve for 7. In all cases, we compute the range for 7 obtained if one 
takes r t /r' t = 1.0 ± Ar, for several values of Ar: 0.01, 0.05, 0.1, 0.2, 0.25. 

The results are shown in Table |l[ (We ignore the discrete ambiguities present 
in the extraction of 7 from Eq. (134]).) For r t /r' t = 1.0 ± 0.25, which is almost the 
full allowed range of r t /r' t , the theoretical error on the extracted value of 7 is about 
20-25%, which is quite large. Thus, as things stand, this method cannot be used to 
make a precision measurement of 7. 

Of course, if the theoretical uncertainty on r t /r' t could be improved, this would 
greatly reduce the error on 7. For example, as shown in Table [L], if the uncertainty 
on r t /r' t were 5%, the error on 7 would only be 2-3°, which would be quite ac- 
ceptable. One possibility for reducing this uncertainty might be to use similar final 
states. Recall that we have assumed that the allowed ranges for r[ had) and r^ ft 
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30 


U.Ul 


27 


o no 

- 32 




a ao 

0.02 


25 


O AO 

- 39 




0.03 


24° 


- 90° 




0.05 
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60 


A A 1 

0.01 


r /i o 

54 


- 64 




A AO 

0.02 


51 


- 74 




0.03 


49° 


- 112° 




0.05 


45° 


- 121° 


80° 


0.01 


73° 


- 84° 




0.02 


69° 


- 95° 




0.03 


66° 


- 125° 




0.05 


62° 


- 130° 



Table 2: The range for 7 calculated from Eq. ([34]) assuming that r t /r' t = 1.0 ± Ar, 
for the input parameters of Eq. (f5§), and for three input values of 7 (ji n = 30°, 60°, 
80°). 

are completely independent. That is, for each of r[ had ' > and , we have assumed 
that the momentum transfer q\ v and the nonfactorizable correction x may take com- 
pletely different values in the two decay processes. However, as noted previously, it 
is quite likely that the momentum transfer q 2 av and the nonfactorizable correction x 
will take comparable values in two similar processes, which will greatly reduce the 
allowed range of r t /r' t . 

Unfortunately, this does not help to reduce the error on 7. Suppose that, instead 
of the values given in Eq. (|35|), the theoretical parameters take the following values: 

Pa = 1.1 , Put = 1.0 , V' ct = 1.15 , V' ut = 1.05 , 

A = 30° , A' = 40° , (3 = 20° . (36) 

Note that the primed and unprimed parameters are similar to one another, so these 
could represent B decays to two similar final states. We again consider three values 
for 7: 30°, 60° and 80°. As before, we can calculate the extracted value of 7 using 
Eq. for r t /r' t = 1.0 ± Ar, with Ar = 0.01, 0.02, 0.03, 0.05. The results are 
shown in Table [|. Regardless of the true value of 7, if the theoretical uncertainty on 
r t /r' t is greater than about 2%, the error on the extracted value of 7 is enormous, 
particularly as regards the upper limit. In fact, even for an uncertainty of Ar = 2%, 
the corresponding error on 7 is at least 15%, which is still large. 

An examination of Eq. ([34] ) reveals why the use of similar final states does 
not help to reduce the error on the extracted value of 7. If the two final states 
are similar, one expects that the experimental observables will also be similar, i.e. 
X ~ X', Y ~ Y', etc. However, in the limit that these sets of observables are equal 
to one another, Eq. (j3~4"D becomes independent of 7, and reduces to the tautology 
1 = 1. Thus, although the error on r t /r' t may be reduced for similar final states, the 
KLY method breaks down if the states are too similar. The net effect is that the 
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error on 7 for similar final states is actually larger than for final states which are 
quite different. 

From this analysis, we conclude that, for similar final states, we really need a 
theoretical uncertainty of Ar < 1% in order to be able to extract 7 with a reasonable 
precision. Unfortunately, with our present knowledge of hadron physics, it does not 
seem possible to definitively establish that Ar < 1% for a particular pair of B 
decays. We therefore conclude that the KLY method works best for pairs of final 
states whose hadronic parameters are quite different. However, in this case 7 can 
only be measured with a precision of ±20-25%, given our current understanding of 
hadronic physics. 

We must emphasize that our calculations have all been done within the frame- 
work of factorization, in which there is still a great deal of hadronic uncertainty. 
However, there is an enormous amount of ongoing work on exclusive hadronic B 
decays. It may well be that, in a couple of years, we will understand hadronic B 
decays well enough to theoretically predict the value of r t /r' t for exclusive states 
with reasonable precision, even for very different final states. If this happens, then 
the KLY method can be used to obtain 7. 

In practice, however, the KLY method will probably be turned around, and will 
be used to learn about hadronic physics. That is, given an independent measurement 
of 7, along with measurements of two b — > d penguin decays, Eq. (|34]) can be used 
to obtain r t /r' t . This information will allow us to test the various models of hadronic 
B decays. 

To summarize, we have re-examined the technique proposed in Ref. || for mea- 
suring (3 (the KLY method). Their original idea is the following: consider the two 
b — > d penguin decays B®(t) — > K°K° and B®(t) — > 4>K S . If one assumes that one 
knows the value of the ratio of two matrix elements in the first process divided by 
the ratio of the corresponding matrix elements in the second process, the CP angle 
(3 can be obtained via time-dependent measurements of these decays. 

In this paper we have pointed out that there is an ambiguity inherent in this 
approach, namely that there are three possible ratios of matrix elements one can use. 
Depending on which assumption one makes, the KLY method can be used to extract 
«, or 7. We have used factorization to show that, in fact, it is the assumption 
which allows 7 to be obtained which is the most accurate. Thus, it seems that the 
KLY method can be used to obtain this CP angle. 

Unfortunately, 7 can not be extracted very precisely. There is a theoretical 
uncertainty in the value of the ratios of matrix elements for the two decays, which 
we estimate to be as large as 25-30%. This leads to an error on 7 of about 20-25%, 
which is substantial. Given our current understanding of hadronic physics, it does 
not appear possible at present to reduce this error. 
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